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BACKGROUND OF THE INVENTION 
1. Field of Ite laventioD 

The presrait innovation relates -to goietally to the methods applied in monitoring of 
physiological fimctions, in particular to methods that are aimed to describe the 
physiological state of the human on the basis of ECG and heart period measurement More 
specfficaUy, the innovation relates to a procedure of d«iving infonnation on the 
physiological state and a measiire of stress of the user on the basis of ambulatory heart 
beatmeasurranent 



Hie pqpsent innovation describes a procedure of segmenting physical state and deriving a 
measure of stress fiom ambulatory heart period signal , 

2. Description of the Prior Art 

Heart period is among the most commonly used parameters in physiological monitoring. 
The wide use of heart period is related, on the one hand, to the availability of 
electrocardiograph (ECG) acquisition device for noninvasive monitoring and, on the other 
han4 to central role of heart period in the autonomic nervous ^'stem fimction and 
sensitivity to several physiological states and conditions. Heart period (or, its reciprocal 
heart rate) forms a basis for different types of analyses and may be defined as the series of 
intervals between consecutive QRS-waveforms in Uie ECG-signal. Another method of 
daiving information on the time distance between consecutive heart beats is the detection 
of heart beat intervals fiom heart pulse signal 

The feet that heart period is a complex product of several physiological mechanisms poses 
a challaige to the use of heart paiod in applied contexts. This is especially the case within 
ambulatory measurement, that is, measurement that is perfonned wifliin natural, fiee- 
Uvmg condition, outside of controUed laboratory environment and protocols. However, the 
muWdeterinined nature of the heart period may also potentate a derivation of additional 
physiological measures fiom the heart period signal by means of decranposing heart 
period mto separate components fliat have a physiological interpretation- 
It is weU known ftat bofli branches of the autonomic navous ^em (ANS), ±c 
sympathetic nervous system (SNS) and parasympathetic nervous system (PNS) influence 
heart rate. It is commonly known that the activity of the SNS and PNS produce 
respectively, an increase and decrease in the heart rate level. It is therefore no surprise that 
much of the work on assessing physiological fimctions and states using infonnation on 
heart beat signal often addresses changes in heart beat as stemming from the influence of 
SNS and/or PNS. Unfijrtunately, it is usuafly very difBcult to determine precisely the 
effects of SNS and PNS on the heart rate, since it is often not apparent which brandi of the 
ANS deteimines changes in heart rate and, in addition to these mechanisms, thee are 
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several other medianisms altmng the level of heart rate hoQi directly and iiulirecfly, matQ^ 
of vMdti are not well-kno\va 

The prior art has documented several research lines atimipting to use heart rate variability 
(HRV) to quantify more selectively the activity of SNS and PNS. It has been documented 
lhat, especially the amplitude of the so-called high fiequency (HP) component of the HRV 
in HoQ fiequency region of 0, 15-0.50 Hz provides information on the level of 
parasympathetic outflow to the heart. Unfortunately, although it has been claimed in some 
instances that the so-called low-fiequency (LF) component of the HRV in the frequency 
region of 0.04-0.15 reflects SNS activity, the effects of SNS on HRV are rather imclear 
and it is known that several ofli^ medianisms also influence HRV and especially Ifae LF 
component, including PNS, hormonal responses, metabolic adjustments, and blood 
pressure control. Thus, increase^ and decreases in the level of heart rate and HRV may be 
due to several sources and therefore, it may be only possible to interpret changes in ihe 
heart rate level and heart rate variability as being indicative of the activity level of SNS 
and PNS during controlled situation and preferably with the aid of other measures. 

The concept of stress refers generally in physiological domain to a state of heightened 
level of physiological activity without immediate or apparent requirements for such 
arousal. In this document, we define a state of stress indicates a body balance wherein tiie 
overall activity level, as indicated in, e.g., cardiovascular output, is substantially hi^er- 
than the level that is required by immediate physical metabolic requirements. The 
physiological state of stress may be due to different sources, such as, for example, 
plQ^sical load (e.g., posture), physical condition (e.g., fewer), mental stress, or emotional 
arousal. 

Feedback and information on personal stress state and more goierally, resources would be 
very helpful for many individuals to monitor and manage their stress levels, to avoid a 
state of burnout, and generally to maintain and enhance health conditioiL Rest and 
relaxation is an important feature of stress management; as it helps to reduce stress and 
fiirth^ buffer against the onset and adverse effects of stress. 

It has been well-documented in the scientific literature that a state of stress is associated 
with heightened SNS influence to the heart and lowered or diminished PNS influence to 
the heart (e.g., Porges 1992). It is also known in the prior art that, at rest during steady 
conditions, relaxation is shown as lowered level of cardiovascular activity and in specific, 
a decrease in the level of heart rate and an increase in the magnitude of the HF component 
of HRV is often found to associate with state of increased relaxation. Some prior work has 
been documented to take advantage of tiie role of HR and HRV in stress and relaxation 
related phenomena (US Patents 4,832,038; 4,862,361; 5,891,044; 5,941,837; 6,104,947; 
6,212,427; 6,358,201). 

Despite this correlational relationship, there has been not veiy much progress in the 
detection of stress-related physiological states on the basis of heart period signal. There 
has been some prior work on using information on heart period and HRV to classify user 
states, in particular in combination with other physiological measures. The prior work 
based on ECG acquisition has been focused mosfly on the detemoJnation of clinical 
condition with special autonomic nervous system tests and is therefore very limited in 
their application to characterize normal behavioral and physiological states in connection 



with, for example, ambulatoiy measurement (US Patents 5,299,199; 5,419,338; 6,390,986; 
6,416,473). 

Th^ has been also some documented work on Hie use of ECG and heart period derived 
measures to detect certain physiological conditions, wherein typically one or more 
parameters are monitored and a threshold limit is set to signal a change in state (US 
Patents 5,267,568; 6,126,595; 6,358,201). These solutions are necessary limited in the 
content of classifying states and suffer fix>m the feet that the signal valiie of ttie heart 
pOTod and HRV parameters is not always the same but mther, typically varies in 
combination with physiological states. In another words, they do not account for the state- 
varying (conditional) relationships between heart period, HRV parameters, and 
physiological states. 

TTiere has been some work ori flie modeling of state- varying relationshq}s. Howevo:, the 
prior work is typically not related to the detomination of stress, may involve heart mte 
measurement but require the use of two or more physiological measures (US Patents 
5,810,0^4; 5,846,206; 5,902,250; 5,921,937). It is thus clear that, fixm the point of 
dififerentiating diffiient user states and in comparison to the acquisition of only one signal, 
these approaches require more eflfort on the measurement of physiological signals and are 
therefore susceptible to involving more mateial costs and more restricted end-user 
protocols. More importantly, the referred work does not include any contribution to the 
identification of stress and relaxation, wherein the occurrence of physical activity has not 
been able to take account in the context of using heart beat signal as a single input 

As indicated above, the major problem in the operationalization, measurement and 
monitoring of stress using information on cardiovascular function, such as acquisition of 
ECG and heart beat, would be the detection and differentiation of the sources of decreased 
and increased cardiac function. This is especially evident with increased cardiac activity 
(e,g., as shown by increased level of heart rate and decreased amplitude of HRV), which 
may result, for example, fix)m increased state of stress, increased state of physical activity, 
or postural changes. 
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OBJECTS AND SUMMARY OF THE INVENTION 

The object of the present invention is provide an automatic analysis tool for the purposes 
of providing information on the physiological state of the user on the basis of ECXj or 
heart period measurement More specii5cally, the object of the invention is to provide a 
procedure for the dififerendation of different physiological states and especially to provide 
a state-detection based method of deriving a measure of stress on the basis of heart period 
data. 

The invention is based on several computational steps wherein the order of flie 
computations has some constraints. The steps may be characterized as follows: (1) initial 
transformations of ECG and/or heart period signals; (2) segmentation of the heart period 
data into stationary epochs; (3) detection of the epochs associated with other-than-stress 
related increases in cardiac activity, including physical exercise, physical activity, 
recovery fiom physical activity, and postural change; (4) detection of segments 
characterized by relaxed state; (5) detection of segments containing a potential stress state; 
and (6)^ combining infomiation obtained in steps 3—5 to provide an overall index of stress. 
The procedure may also contain an initial set-up of parameters wherein some properties of 
the state detection system are either inputted or determined automatically. 

According to the present system, the detection of stress , state is based on tiie following 
physiological assumptions: There is sympaftietic dominance in relationship to the 
parasympathetic dominance and there is no evidence of physical activity, exercise, 
movmioit or posture influmced cardiac reactivity. 

The invention may be described as an expert system that consists of a sequence of 
computations and inferences and provides a novel approach to the detection of 
physiological states and especially stress on the basis of heart beat signal. The types and 
methods of giving user feedback are dependent on the purpose of the equipment and 
sofiware in connection with v^ch the classification and state detection procedure has 
been applied. 

The invention may be applied to and in association with devices such as heart rate 
monitors and otiier wearable and mobile computing devices, other types of equipments for 
physical monitoring and especially involving ECG detection, and software products suited 
for the analysis of heart p&nod signal. The present invacition is usefiil in combination of 
any physiological monitoring in flie area of monitoring, enhancing and optimizing 
physiological resources to better manage healthy lifestyle and wellbeing, sports training 
and fitness, and working capacity. It may be especially usefiil in the monitoring and 
providing information on state of stress, vsdi^ it may be used in the long-tmn mcmitoring 
of the resources. 



BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 . An overall view of the sequences in the procedure. 
Figure 2. Initial trans&rmations of the ECG and heart period signal 
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Figure 3. An illustration of the procedure for the detection of exercise, recovery from 
exercise, physical activity, movOTient and postural dianges. 

Figure 4. Procedure for the detection of relaxation, stress state, and combination of 
information to desoibe overall body resources. 



DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

INTRODUCTION 

The innovation is described here with the aid of an example implementation. It should be 
noted that the described system is not bound to any specific model or specifications, but 
rathor, diJSerent alterations, fomis, and improvements are possible and are in line with the 
spirit of the innovation. Thus, Uxo following ma:ely contains a desoiption of flie preferred 
embodiments of the innovatioiL . 

The system invoh^es initially an estimation, input, or use of previous values to characterize 
individual parameter values that influence the detection of states within the system. This 
estimation of the values is optional for the system. The values may be estimated by using 
fonnulas based on empirical data and user-iriputted values, such as age, weight, height, 
and sex. The values may be also based on history values based on empirical physiological 
measurements, or they may be iiqjutted by the user or expert. Examples of values that may 
be used on this context are maximum and minimum levels of OTQrgen consumption, heart 
rate level, or HRV level. 

Figure 2 shows an overview of the initial transformations of the ECG and heart period 
signal. The ECG-signal is transformed into consecutive R-R intervals by using sequential 
QRS-complexes as markers of the beginning and end of a R-R interval. Also other 
components of the ECG signal may be used to detect the RR interval. The consecutive R 
R intervals are scanned through an arti&ct detection filter to perform an initial correction 
of felsely detected, missed, and premature heart beats. Another potential methods of 
deriving information on the intervals between consecutive heart beats would be based on 
the detection of intervals fiom pulse signal on the basis of different detection algorithms. 
It should be therefore noted, that in principle, the present innovation applies to all domains 
wherein temporal information is obtained on the timings between consecutive heart beats. 

Consecutive artifact- corrected R-R intervals are transformed into equidistant heart period 
time series by using a weighted linear interpolation of the R-R intervals. For example, a 
time domain sampling rate of 5 Hz may be used. 

Initial features that may be used in the segmaitation of the physiological process of the 
user into stationary epochs are computed prior to the segmentation. The features may 
include several different components, including the following: spectral power (ms^) in the 
LF and HF firequency regions, information on respiratory period and ventilation as 
derived, for example, from flie heart period according to Kettunen & Saalasti (Patent 
pending, FIN 20025029), variation and inconsistencies in respiratory parameters, and 
oxygen consumption as derived from, for example, heart period and heart period derived 
respiratory period or ventilation. If available, also external measures of physical activity 
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may be used as a component for segmenting Ihe physiological processes into stationary 
epochs. All these measures are computed on a continuous basis with preset or dynamically 
alternating window sizes, thus potentiating the capabiUly for lealrtime use. 

A univariate or multivariate set of features is combined to segment the physiological 
processes into stationary epochs. In another words, the target of the present procedure is to 
divide segments wherein Ihe properties of Ihe physiological ^stmi are similar and 
consistent witihin-the-segment, using the information obtained from the derived 
physiological features. An example of a algorithm to perform this is the so-called 
generalized likelihood ratio test, which basically poforms the minimization of the 
variance wittiin ttie s^^ients and maximization of the variance between the (consecutive) 
s^ments. This process is controlled by threshold parameters that determine the sensitivity 
of Ihe segmentation processl fhe use of multivariate set of parameters describing user 
state is recommended to gain more stability and reliability to the segmentation of data into 
stationary epochs. 

The g^eralized likelihood ratio test is classically applied to model selection but it is also 
adapted to segmentation algorithms, e.g. the one presented by Fancourt and Principe. In 
this embodiment the fiinction used to estimate the signal in the segmmt is median. The 
error is calculated fiom the mean absolute error between the median and tiie sigtml within 
segment For clos^ details of the algoritiun see the article by Fancourt et al. 

In Ihe following the procedure for detecting physiological state is presented Selected 
physiological parameters are analyzed from each stationary segment to form a basis for 
inference on states that include increased cardiac activity. The detection of physically (i.e., 
metabolism related) induced states of increased cardiac reactivity is excluded before the 
detection of potmtial states and intCTsities of relaxation and stress. Figure 3 presents a 
flowdiart on the procedure of detecting physical activity related states. 

Increased level of cardiac activity (e.g., increased heart rate level or decreased heart period 
level) may be due to frie process of recovering from physical exercise or from any physical 
activity. The state of recovering from exercise may be detected by, for example, applying 
the method described in Saalasti, Kettunen, PuDddnen (A procedure for predicting body 
fetigue during exercise and recovery from exercise on the basis of physiological 
measurement. Patent pending, FIN 2002) to deteiinine the level of recovery requirements 
in the body arid classify seginOTts bvolving recovery demands higher than a cotain 
predetermined threshold as recovery state. 

The syst^ requires information on the intensity of physical activity and exercise to 
differentiate raercise related effects on cardiac system from, e.g., stress related effects. 
The transformation of heart beat level into proportional intensity of physical activity, 
applying for example information on the relationship of heart rate level to oxygen 
consumption. It £ possible also to use information on initial parameters, such as maximum 
and minimum levels of heart rate and oxygen consumption to increase the accuracy of 
detamining the intensity of physical activity. However, it is clear that the direct 
transformation of heart level to exercise intensity does not produce optimal results for the 
purposes of differentiating other tiian exercise related increase in cardiac activity fiom 
exercise induced reactivity. 
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Tlie estimate of the proportional intCTsity of physical activity may be enhanced by 
including infonnation on the ventilation, as obtained, e.g., from heart period. Furfliermore, 
given that ex^cise has also certain pattern of changes in respiratory activity (e.g., exercise 
is associated with increased respiration, whereas otfaer-than-exercise related cardiac 
reactivity is not often associated with similar respiratory reactivity) and HRV (e.g., a 
certain level of HRV is typically associated with certain level of heart rate during 
exmnse), and temporal length (i.e., exercise and intensive physical activity is necessarily 
bound to have a certain temporal length). It woxild be also possible to use external sensors 
to detect the occurrence and intensity of movement and physical activity with, for 
example, accelerometer or skin temperature sensors, to support the determination of 
exercise int^osily. 

A certain threshold value is ^preset to determine whether a certain segment involves 
exercise or intensive physical exerdse. Given that no ottier reasons related to respiration, 
HRV, or temporal length of increased heart rate provides evidence that a potmtial exercise 
induced increase in cardiac activity is not associated with intensive physical activity, 
segmeiijs containing a mean intensity of {diysical activity hi^^ than, for example, 50% 
of exercise intensity may be classified as int^ive exercise. 

Movements, start-up of a physical activity, and postural dbanges have all metabolic 
requirements and also ina*ease heart rate level. It is thus necessary to differentiate tiiese 
effects from other, non-metabolic factors influencing cardiac activity. It is known that 
PNS often controls increase in cardiac activity at the onset of movement and decrease in 
cardiac activity at the end of movetnent It is also known that standing up, which poses a 
so-called orthostatic reaction, is associated with increase in the spectral power in the LF 
fiequency region of HRV in relationship to HF power, which usually decreases when 
standing up. Based on these assumptions, a covariance parameter according to Equation 1 
has beai datived to index movement related changes in cardiac activity. 

Equation 1 



movement = oc* cov{}og{HFpow)) + 




An increase in the level of the described parameter indicates the prevalence of mov^ent 
related inoieases in cardiac activity and thus may be used to differentiate states of 
movement, onset and offset of physical activity, and postural changes. 

The detection of physical activity is similar to that of intensive physical activity but 
the threshold values are different However, the detection of light physical activity differs 
from that of detecting intensive physical activity in that it is expected that metabolic 
demands related changes in cardiac activity are associated with light physical activity. In 
another words, the covariance t^m introduced in Equation 1 is used to differentiate light 
physical activity related increases in metabolism from cardiac reactivity that is due to 
differ^ sources. It is also h^ of note that, in a similar manner to the int^isive physical 
activity, additional sensors detection movanent could be also used to differentiate 
segments that include physical activity. 

The detection of movement and reactivity due to changes in posture is performed by using 
the covariance term presented in Equation 1 with preset threshold values and a period 
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sunounding the covariance. The covaiiance term be smoothed with, for example, a 
Hanning window. 

In general, it is known that exercise, physical activity, and movement all have their typical 
time fiames that have minimal temporal requiiemaits for recovery. For this purposes, it 
may be possible to vise a frequency or time domain measure to provide information on the 
temporal properties of the changes in cardiac activity, thus potentiating the comparison of 
those with those ^ical in the case of differrait types of physical activity. 

It is important to notice that, as for example, the changes on the heart beat during physical 
activity are associated with seyeral differrat types of changes, ftese changes may be 
combined in an automatic decision fimction that may be, for example, detaministic, 
heuristic^ or based on multivariate methods such as fiizzy systems and neural network. 

Figure 4 shows an overall view on the detection of relaxation and stress. Relaxation index 
for the segment is determine by the combination of heart period and HF power as 
illustrated in Equation 2. 

Equation 2 



, ,j iiR max- /fl? min . rrry- 
relax _ threshold = — + JnUi mm 

t^seg men t end 

relax _index {segment) —irffsmertt^^^fri 



HFpowjt) 



HR{t) 



segment _ length 



The relaxation may be detected for example by calculating mean heart rate in the segment 
and flien using a threshold proportional to minimum and maximum heart rate. An example 
is presented in Equation 2. If mean heart rate is less than the threshold the segment may be 
detected as relaxation. 

If no relaxation state is detected, the data segment is a potential candidate of including 
stress. Equation 3 shows an example of how an index of stress may be computed. In 
Equation 3, HR denotes heart rate level, CT denotes inconsistencies in the frequency 
distribution of HRV due to changes in respiratory period, or altematively, variability in the 
respiratory signal HFpow and LFpow denote spectral powers in the HF and LF regions of 
the HRV, respectively. 

Equatton 3 

t^sesm^^end HR{t)'CT{t) 

. ^ , t^eJ;i s.aMHFpOWit))'\O^LFpOMit)) 

stress index(segment) = = ; — 

~ segment length 



Of course, both indices shown in Equations 2 and 3 are only illustrative and different 
fomis may freely apply. Both indices may be also used to derive continuous parameters. 
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vMda may be useful in some contexts of deriving ovei:alI, long-time infonnation on stress 
and relaxation. 

It should be clear to anyone femiliar with fte field that the formula presented for 
computing stress index is only illustrative and may be formed in several methods. The 
principle is to search for periods of increased heart rate witii other markers of potential 
stress and related in general to heightened level of SNS activity and diminished level of 
PNS activity. 

The indices of relaxation and stress may be combined taking advantage of tiie proportion 
of diffeient states and the intensity of relaxation and stress states. An example of sudi 
combination is presented in Equation 4. 

Equation 4 

total = (1 - RLXpow^ 100) + C '^'STItpow, 

where c is a constant, T is total time of the measurement, Tj^is time classified as 
relaxation, is time classified as stress, RLXpow the intensity of relaxation in 

percentages and STRpow is the intensity of stress state (see the computation of stress index 
in Equation 3). The inteaosities may be calculated via mean or median values of the 
corresponding time series indices. The combined index may be especially usefid in flie 
comparison of different days. 

In principle, a segment is identified as stress if heart rate is constantiy at a high level and it 
is apparent that it is caused by lowered level of PNS, higher level of SNS, and there is no 
evidence of physically detamined heart rate reactivity. In feet, given tiiat stress is 
associated with a tonic low level of PNS and that high level of PNS potentiates reactivity, 
the absence of PNS mediated responses is also an indicator of stress and may be taken 
more formally into account. If a person is highly stressed, it may not be possible to find 
states of relaxation even during long pmods. 

Feedback on the results of the present procedure may be presented to tiie user either on 
real-time ^plication or off-line, after the measurement. It is possible to give feedback on 
states, on relaxation and stress components, and overall resources. The feedback may be 
giv^ in many forms, for example, in graphics or as a single parameter that is easy to 
interpret 




Claims 

What is claimed is: 

1. A procedure of segmenting heart beat signal into physiological states, 
characterized by at . least one of the following 

a. A method for the segmentation of the heart beat signal mto intemally 
coherent segments 

b- The use of at least one analysis method for the identification of segments 
that include mcreased metabolic rate due to exercise, physical activity, 
movement, or postural changes, including methods such as, for example, 

i. Information on the co occurrence of heart beat changes and HRV 
measure?, such as the moving covariance presented in the preferred 
embodiments 

ii. Information on the HRV and any component of it, such as, for 
example, LF and HF components 

iii. Exercise intensity, such as that based on tiie transformation of heart 
beat level and/or information on respiratory period into level of 
oxygen consumption 

iv. Recovery from exercise 

V. Respiratory period or ventilation in combmation with heart rate 
level and/or HRV or its decomposition 

vi. The use of the information on the temporal properties of the 
exercise, physical activity, and movement as characterized by, for 
example, frequency or time domain methods 

vii. The use and combination of information as obtained from several 
detection 

2. A procedure according to Claim 1, which is used to identify states and periods of 
exercise, physical activity, movement, or postural changes from heart period signal 

3. A procedure for the detection of stress and resoiirces, characterized by a procedure 
according to Claim 1 and combining the procedure of Claim 1 to the information 
derived on the basis of heart beat data, based on, for example, 

a. Heart beat level 

b. Derivation of index or measures of stress based on information obtained 
from the heart beat 

c. Derivation of index or measures of relaxation based on information 
obtained from the heart beat 

d. Derivation of index or measxires of resources based on information 
obtained from the heart beat 

4. A procedure for the computation of overall index of stress, relaxation and 
resources, characterized by the procedure according to Claim 3 and by the 
computation of an index that represents a suirmiary of the existence and level of 
stress, relaxation and resources for a period of measurement, such as a 3 hours, 
half day or ftdl day 
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5. A system for measuring stress and relaxation on the basis of heart period 
measurement, wherein information on the length of detected relaxation and length 
of detected stress is used as informative in the detection and quantification of 
relaxation and stress states 

6. A system for measxiring stress and/or relaxation on the basis of heart period 
measurement, wherein the measure of stress and/or relaxation is computed for a 
segmented period of heart period 

7. A system according to the Claims 1 5, wherein information on the exercise level, 
physical activity, movernient, or postural changes is obtained from heart beat and/or 
one or more sensors detecting movement, such as accelerometer 

8. A system according to claims 1 3, which is used for correcting heart rate based 
oxygen and energy consumption estimate 
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ABSTRACT 

The present innovation is a procedure for the segmentation and classification of the user's 
state and for the determination of a measure of stress on the basis of ambulatory 
measurement of heart beat. The procedure may be described as an expert system that is 
based on several computational steps. The system provides information for tiie automatic 
detection and differentiation of movement from other sources of cardiac reactivity and 
performs a detection of states with the capability of functioning in ambulatory conditions 
wherein heart beat measurement often contains complex components. The core of the 
innovation is the differentiation- of the source of increased cardiac activity between the 
stress and physical activity induced accelerative changes in the cardiac function. The 
information obtained from the classification of physiological state may be used as a 
starting point in analysis of aiyibulatory ECG or heart beat information and directly to 
provide the user simple feedback on the physiological state and especially on a 
physiological state of stress. The procedure has wide spread use in the ambulatory 
monitoring of physiological signals and states, and especially as related to the devices 
aimed at providing feedback on users health and wellbeing. 
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